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Abstract 



We proposed a model of interacting market agents based on the 
generalized Ising spin model. The agents can take three actions: 
"buy," "sell," or "stay inactive." We defined a price evolution in terms 
of the system magnetization. The model reproduces main stylized 
facts of real markets such as: fat-tailed distribution of returns and 
volatility clustering. 
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1 Introduction 



Financial market modelling lies in the field of interests of both theoreticians 
and practitioners. Among those who develop these models are also econo- 
physicisits. A view on a financial market as a complex system of investors 
similar to complex physical systems appeared to be very fruitful, since it al- 
lowed to reproduce many characteristic features of such a market [H El EJ H] . 
A special group constitute approaches based on the Ising model [HJ El 
[?l [8j or its generalization to a three-state model P, HO, HT] . They identify 
agents with spin variables which can take specific values depending on agents' 
decisions. In much the same way as Ising spins, agents interact with each 
other, which leads to a herding behavior and, as a consequence, bubbles or 
crashes. An important target for such models is a reproduction of real market 
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stylized facts such as fat-tailed returns distribution, clustered volatility, or 
long range correlations of absolute returns |12j . 

In this paper we proposed a simple model of financial markets, based on 
the Granovetter threshold model of collective behavior, which is in a good 
agreement with the above facts. 

2 Description of the model 

An inspiration for our work was the Bornholdt model [13], corresponding 
to the Ising spin model with an additional minority term, where agents act 
under the influence of their neighbors (Ising part) and a global magnetization. 
The local field for i-th spin in the Bornholdt model is defined as: 

N 1 N 

hi^^JijSj-as^ — ^s^, (1) 

3=1 3=1 

with a global constant a > 0. The model interprets the price of an asset 
in terms of the magnetization, which enables the authors to reproduce some 
stylized facts of financial markets. 

Here we developed a generalization of the Ising spin model that also uses 
the absolute value of the magnetization as a factor controlling the dynamics. 

Let us consider a model of N interacting market agents where each agent 
takes one of three actions: "buy," "sell," or "stay inactive." Such an agent 
can be represented by a three-state spin variable st(t) taking values +1 when 
the agent is buying, —1 when the agent is selling, and in the remaining 
case. 

The agents interact with each other according to an interaction matrix 
Jij. In our model we assumed the interaction matrix corresponding to a 2- 
dimensional squared lattice in which each agent interacts only with its four 
nearest neighbors with an equal strength J. The interaction strength can 
be ferromagnetic (J > 0) when investors try to act like their neighbors or 
antiferromagnetic (J < 0) when they try to play against their neighbors. We 
chose the ferromagnetic case to introduce the herding behavior. 

At each time step t the agent % takes its value according to the following 
formula, 

N 
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where sign is a threshold signum function, 

1 if x > q, 
sign q (x) = ^ if — q < x < q, (3) 

— 1 if x < —q. 

The function r]i(t) has random values from the Gaussian distribution with 
mean and variance equal to 1. The term arji(t) simulates an individual 
erratic opinion of the i-th investor and the parameter a defines strength of 
individual opinions. 

We define the magnetization of the network, 

1 - 

M W = ^5>W> ( 4 ) 

i=l 

which together with a constant A forms a threshold parameter q of the sign 
function. Let us notice that for A = our model is identical with the 2- valued 
Ising spin model. 

We simplified the procedure of price calculation presented in [T3] by omit- 
ting the influence of chartists and fundamentalists in the population of in- 
vestors and redefining the price of an asset as: 

P(t) = P (t)e M(<) , (5) 

where Po{t) is, in agreement with the efficient market hypothesis |14| . a 
geometric Brownian walk corresponding to fundamental price changes. Let us 
put Po{t) = Pq constant for more simplicity. We got therefore a logarithmic 
rate of return: 

r(t) = M(t)-M(t-l). (6) 

According to (2), each agent is under the influence of three factors. The 
first is an imitation of their neighbors, associated with the matrix Jy which 
is responsible for the herding behavior. The second factor is an individual 
opinion of the agent provided by the term ar)i(t). So far, it is the standard 
Ising model. Yet, there is one more factor, X\M(t — 1)|, which plays the 
role of a threshold parameter. Only those agents that are able to exceed 
the threshold are allowed to trade. The value of the threshold depends on 
the absolute magnetization. According to (5), the magnetization measures 
a deviation from the fundamental value. So, when it is large, the agents are 
afraid of trading, unless they have a strong support from the neighbors or 
from their private opinions. 
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Figure 1: Returns r(t) in time t for parameters J = 1, a = 1 and different A. 



3 Results of the simulation 



Using a square 32 x 32 lattice of agents with periodic boundary conditions 
we computed the history of the magnetization, and thus the evolution of 
the stock price. The agents were allowed to interact only with their nearest 
neighbors. 

The simulation was started with a random configuration of spins. At each 
time step a randomly chosen spin was updated according to the evolution 
equation (2) and this was repeated N times. The first 5000 time steps were 
ignored as a period of system thermalization. 

We observed that for proper parameters price returns show volatility clus- 
tering (fig. 1) and fat tails (fig. 2). 

For real prices, fat tails of returns distribution (r r (t) = \og(P(t)/P(t— r))) 
are getting slimmer with rising time lag r. Finally, the distribution became 
normal for sufficiently large r. In Figure 3 we presented the distribution r T (t) 
for different r. The tails are slimming down with rising r. The distribution 
has Gaussian tails although its shape is slightly different. 

Volatility clustering can be quantitatively shown using the autocorrelation 
function of a time series x(t): 



'-. i: 
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(7) 



where a x is the variance of x(t), and (...) means the average over t. We com- 
puted the autocorrelation function of returns C r (t) (fig. 4) and the absolute 




Figure 2: (Color online) Histogram of returns rie(t) for parameters J = 1, 
o = 1 and different A (left picture), and for parameters J = 1, A = 10, o = 2, 
and cr = 0.5 (right picture). 
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Figure 3: (Color online) Histogram of returns r T (t) for parameters J — 1, 
a = 1, A = 10, and different time lag r. The returns were rescalled by the 
corresponding standard deviation a r . The distributions were shifted upward 
with the factor 10, 100, 1000, ... The normal distribution was plotted with a 
dashed line. 
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Figure 4: Autocorrelation function of return C r (r), obtained for parameters 
J = l, a = 1, A = 15. 



value of returns C\ r \(t) (fig. 5). The autocorrelation of returns decays very 
fast, which is consistent with observations of real markets. It has also the 
same shape as the distribution presented in [10]. The autocorrelation func- 
tion of the absolute returns is a slowly decaying exponential similar to the 
results of [13] and [10]. 



4 Conclusions 

We proposed a simple model of financial markets based on the Granovetter 
model. We introduced a threshold controlled by a magnetization of the 
system. This makes an agent take an action only if its confidence is strong 
enough to overwhelm the threshold. 

We defined a logarithmic rate of return as a change of the magnetiza- 
tion. The model with such a defined price reproduces main stylized facts of 
financial markets meaning a fat-tailed distribution of returns, volatility clus- 
tering, very fast decaying autocorrelation of returns and much slower decay 
of autocorrelation of absolute returns. 

It has been observed for real markets that the distribution of returns 
r T (t) becomes Gaussian for large r. The distribution of returns generated 
by our model approaches the Gaussian distribution for rising r, however its 
shape differs from the normal distribution. The autocorrelation function 
of the absolute returns is a slowly decaying exponential function, while the 
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Figure 5: Autocorrelation function of absolute return C\ r \(t), obtained for 
parameters J=l,er=l,A = 15in log- log scale, and in semi-log scale 
(inset). 



empirical study of the markets reveals a power-law behavior. Both these 
issues can be considered as weaker points of our model. 
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